} t 


Department  of 
Aerospace  Engineering 
and  Applied  Mechanics 


ft 


SURFACE  RADIATION  INTEGRALS  IN  AERODYNAMIC  SOUND 


Final  Report 


by 


MORRIS  P.  ISOM 


November  1974 


U S.  ARMY  RESEARCH  OFFICE 
GRANT  No.  DAHC  04-74  G-0077 

D D C 


Approved  For  Public  Release. 
Distribution  Unlimited 


POLY-AF/AM  Report  \o.  74-PO 


The  findings  in  this  report  are  not  to  bo  construod 
as  an  official  dopartmont  of  tho  army  position,  unloss 
so  dosignatod  by  othor  authorizod  documents 


r 

I 


line  1 a r»r.  it'jr  i 

StCUntV  Classification  _ 

DOCUMENT  CONTROL  DATA  - R & D 

{Security  cle*»iltt  etron  ot  title,  body  ol  abetrmci  and  inJentng  annotation  thuet  be  ottered  when  ih^ovorel^ogot^e^^eeoUlod^ 
ORIGINATING  ACTIVITY  (Coipo.Al.  ai»a>)  jla.  «»0«T  IlCum'*  CGAttlFICATION 

Unclassified 

Polytechnic  Institute  of  New  York  group  “ 


Surf, ice  Radiation  Integrals  in  Aerodynamic  Sound/ 


• tifU  «iPf  vE  NOTH|  Tyne  ol  report  and  Incluolve  date 

Final  Report  1 ' ’ / I / 7 3 - 5/3'/74 


991 


fa.  TOTAL  mo  OF  PA6KI 


1b  NO  OF  RIM 


d\S-k-£/Ai\\-n~'Jf 


Oa.  OAlOlKATO***  REPORT  NuMRIRHj 


noly_AF/AM  Report  No.  74-20 


TH|R  RtPOR  r NO(l)  (Any  other  number  a that  mar  be  i 
' a report) 


11  IPONSORINC  MILITARY  ACTIVITY 

U.  V, . Army  Research  Office 
Box  CM,  Duke  Station 
Durham,  North  Carolina  27706 


Equations  are  derived  for  the  contribution  to  far  field  acoustic  pressure  from 
sources  distributed  on  a rotating  and  translating  rigid  surface,  and  from  shock  waves 
carried  by  the  surface.  Radiation  integrals  are  expressed  in  coordinates  fixed 
relative  to  the  accelerating  surface.  The  main  results  are  an  extension  of  the 
generalized  surface-fixed  Lowson  formulation  of  the  aerodynamic  sound  problem  as 
derived  by  Ffowcs  William  and  Hawkinqs  in  1968.  The  qeneral  formulation  is  intended 
for  application  to  the  problem  of  calculating  the  sound  field  from  a helicopter 
rotor  blade  with  transonic  conditions  in  the  blade  tip  region.  The  sound  field 
therefore  includes  contributions  from  non-compact  source  distributions  on  the 
blade  surface  and  shock  waves. 


ABSTRACT 


Equations  are  derived  for  the  contribution  to  far  field  acoustic 
pressure  from  sources  distributed  on  a rotating  and  translating  rigid 
surface,  and  from  shock  waves  carried  by  the  surface.  Radiation  integrals 
are  expressed  in  coordinates  fixed  relative  to  the  accelerating  surface. 
The  main  results  are  an  extension  of  the  generalized  surface-fixed  Lowson 
formulation  of  the  aerodynamic  sound  problem  as  derived  by  Ffowcs  William 
and  Hawkinqs  in  1968.  The  general  formulation  is  intended  for  application 
to  the  problem  of  calculating  the  sound  field  from  a helicopter  rotor 
blade  with  transonic  conditions  in  the  blade  tip  region.  The  sound  field 
therefore  includes  contributions  from  non-compact  source  distributions  on 
the  blade  surface  and  shock  waves. 


IVP  "-I* 

I; it  imte' 

USE  ld<l  SKtlM 

JBSTlf iCIUM  


It  

liSTtilt: 


TABLE  OF  CONTENTS 


Page  No. 


ABSTRACT 

1.  INTRODUCTION  I 

2.  SURFACE  INTEGRALS  AND  AN  INTEGRAL  THEOREM  3 

3.  SIMPLE  IDENTITIES  7 

4.  GEODESIC  PROPERTIES  OF  THE  ROTATING  SURFACE  12 

5.  SHOCKS 

1.  The  hovering  rotor  14 

2.  Advancing  rotor  17 

6.  SUMMARY  AND  CONCLUSIONS  19 

7.  REFERENCES  20 

NOTE  ADDED  IN  PROOF  21 


SURFACE  RADIATION  INTEGRALS  IN  AERODYNAMIC  SOUND 


I.  INTRODUCTION 

The  sound  field  produced  by  a rigid  surface  in  arbitrary  motion 
was  described  in  general  terms  by  Ffowcs  Williams  and  Hawkings  in  1968 
(Reference  l).  Several  alternative  descriptions  of  the  radiative  component 
of  the  far  pressure  field  were  given.  Each  formula  is  an  asymptotic  solu- 
tion to  an  inhomogeneous  wave  eguation  that  is  a generalization  of  the  wave 
eguation  used  by  Liqhthill  to  describe  aerodynamic  sound. 

One  of  the  formulations  of  Reference  1 was  singled  out  by  Farassat 
(References  2 and  3)  as  suitable  for  detailed  noise  calculations.  His 
results  are  a considerable  clarification  of  the  original  Ffowcs  Williams 
and  Hawkings  analysis,  and  he  has  extended  them  to  the  point  where,  in 
pr i nc i pie, they  are  all  that  are  necessary  for  numerical  calculation  of 
a sound  field  that  includes  non-compact  source  distributions. 

The  unusual  and  useful  feature  of  Farassat's  work  is  the  coordinate 
space  in  which  radiation  integrals  are  expressed.  This  space  is  a tem- 
poral and  spatial  subspace  of  the  complete  four-dimensional  space  and 
time  coordinate  system.  The  possibly  singular  Doppler  factor  is  suppressed 
in  all  radiation  integrals  in  this  coordinate  system,  and  this  result 
is  an  important  advantage  of  the  method.  The  elimination  of  an  aerodynamic 
(Doppler)  singularity  from  radiation  integrals  is  accomplished  at  the 
expense  of  a geometric  singularity  that  appears  in  all  radiation  integrands, 
and  this  new  singularity  occurs  at  all  Mach  numbers.  The  singularity  appears 
when  the  retarded  unit  radiation  vector  from  a point  on  the  moving  surface 
to  the  distant  observer  is  parallel  to  the  surface  normal  vector.  Farassat 
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has  shown  in  Reference  2 that  the  singularity  is  integrable  for  a useful 
class  of  surface  geometries.  It  seems  likely  that  it  can  be  accomodated 
in  numerical  calculat ions , but  whether  it  will  be  a minor  or  major  an- 
noyance in  specific  applications  is  not  yet  clear.  In  any  case,  the 
methods  and  results  of  References  2 and  3 may  be  considered  as  alternatives 
to  Lowson%  formulat ion  of  the  aerodynamic  sound  problem  (Reference  4)  and 
to  generalizations  of  Lowson's  work  that  appear  in  Reference  1. 

Here  we  intend  to  return  to  the  Lowson  description  of  the  sound 
field  in  its  generalized  form  stated  by  Ffowcs  Williams  and  Hawkings. 

The  motivation  is  based  partly  on  the  fact  that  the  results  of  References 
1 and  2 strongly  suggest  that  the  Doppler  factor,  which  is  important  in 
low  speed  problems,  appears  to  play  no  role  in  a high  speed  problem.  It 
would  therefore  seem  likely  that  it  can  be  eliminated  from  the  surface- 
fixed  coordinate  description  of  the  sound  probiem  as  originally  given  by 
Lowson  and  generalized  in  Reference  1.  This  conjecture  is  true  to  an 
interesting  and  significant  extent  for  the  problem  of  sound  radiated  by 
a rotating  and  translating  surface.  However,  our  results  have  one  possible 
disadvantage:  while  all  radiation  integrals  have  neither  a Doppler  or 

geometric  singularity  in  them,  they  all  appear  with  time  derivatives  in 
front  of  them.  A considerable  and  successful  effort  was  expended  in 
References  1 and  2 to  get  all  derivatives  of  radiation  Integrals  under 
the  integral  sign.  The  formal  ism  converts  derivatives  with  respect  to 
observer  coordinates  and  time  to  derivatives  with  respect  to  coordinates 
defined  on  the  moving  surface.  These  transformations  are  essential  if  the 
main  purpose  is  to  obtain  rough  estimates  of  radiation  integrals.  However, 
they  have  the  effect  of  increasing  the  strength  of  geometric  singularities 
that  must  be  integrated.  On  the  other  hand,  the  appearance  of  time 
derivatives  before  integrals  with  well-behaved  integrands  does  not  appear 
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to  pose  any  problem  for  numerical  calculations,  and  this  form  may 
in  fact  make  numerical  calculations  easier. 

General  results  will  be  qiven  here  for  a particular  class  of 
radiation  integrals.  These  are  surface  radiation  integrals,  or  integrals 
of  shell  distributions  of  sources  as  they  are  called  in  Reference  1.  They 
make  an  important  and  possibly  dominant  contribution  to  far  field  sound. 

The  general  formulation  is  intended  for  application  to  the  problem  of 
calculating  the  sound  field  from  a helicopter  rotor  blade  with  transonic 
conditions  in  the  blade  tip  region.  The  problem  therefore  includes  non- 
compact source  distributions  on  the  blade  surface  and  radiation  from 
shock  waves  carried  by  the  blade. 


2.  SURFACE  INTEGRALS  AND  AN  INTEGRAL  THEOREM 
Let  S be  the  surface  of  a rigid  body  that  rotates  with  constant 
angular  velocity  u)  about  an  axis  fixed  in  space  and  that  translates 
with  forward  velocity  v in  a plane  perpendicular  to  cu  . The  contribu- 
tion to  far  field  pressure  p from  a surface  distribution  of  sources 
of  strength  F on  S is 


4tt  p 


I*  F cos  ft 
*t  J cr 


(2.1) 


Notation  Is  the  same  as  in  Reference  1;  */^t  is  time  differentiation 
in  a reference  frame  fixed  relative  to  undisturbed  air  and  c is  the 
constant  sound  speed  of  the  undisturbed  medium;  0 is  the  angle  between 
the  unit  radiation  vector,  directed  from  a point  on  S to  the  observer, 
and  the  outward  surface  normal  vector;  r is  the  distance  from  an 
Integration  point  on  S to  the  observer,  and  Mr  is  the  component  of  the 
body  surface  Mach  number  in  the  radiation  direction.  Integration  is 
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over  the  body  surface  S with  the  integrand  evaluated  at  retarded  time 
t = t - r/c . 

Let  7)  = (T]j  »T|g»71j ) be  a Cartesian  system  of  coordinates  fixed 
relative  to  S and  let  x = (x^.x^.x^)  be  observer  position  in  a frame 
fixed  relative  to  undisturbed  air.  Then  F,  cos  9 and  M^  are  generally 
explicit  functions  of  position  7),  x and  t through  their  dependence  on 
retarded  time  t,  which  is  a function  of  T,  x and  t.  To  emphasize 
this  dependence,  (2.1)  is  written  as 


4tt  p = ^ Jf(t  ffl.x.t)  ,fi) 


cos  ft  dS 
cr  (1-M^)  * 


(2.2) 


with  similar  t and  T)  arguments  for  cos  9 and  M^ . 

The  general  representation  of  sound  from  a surface  integral  has 
a Doppler  factor  |l-Mrl  ',  and  thus  includes  radiation  Mach  numbers 
Mr  that  may  be  supersonic.  Here  we  assume  that  < 1 and  exclude 
supersonic  surface  motion.  Of  course, there  may  still  be  regions  of 
supersonic  flow  over  S and  embedded  shock  waves. 

The  integral  in  (2.2)  will  be  transformed  in  the  surface  fixed 
frame  so  that  the  factor  ( 1 — ) * is  suppressed.  The  transformat  ion 

is  based  on  some  identities  and  an  elementary  theorem  of  differential 
geometry.  We  must  first  define  the  operations  of  total  surface  divergence 
and  partial  surface  divergence  relative  to  the  surface  S . 

Let  F be  an  arbitrary  vector  function  of  retarded  time 
t (T , x , t ) and  position  ^ on  S . Let  n(^)  be  the  unit  outward  normal 
vector  on  S specified  in  V,  coordinates,  so  that  n is  independent 
of  t . The  radiation  vector  r(r)  is  the  unit  vector  from  a point  on 
S to  the  observer.  Two  unit  vectors  P(t,1|)  and  m(T,?i),  both 
tangent  to  S , are  defined  as  in  Reference  I by 
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, m = n x f , |?xn|  = sinO  > 0 , 


so  that  T,  m and  n form  a right-handed  triplet  of  unit  vectors  on 
S , specified  in  coordinates  fixed  to  S . The  vectors  f and  in 
play  a fundamental  role  in  Farassat's  work.  Here  they  may  be  considered 
as  vectors  derived  from  the  more  fundamental  unit  vectors  r and  n and 
always  expressible  in  terms  of  ? and  n. 


We  define  the  tota 1 s ur face  d ? verge nee  of  any  vector  F defined 


on  S by 


V,.  • F = f " VF  • f + m • 7F  * m 


The  gradient  operator  V in  (2.4)  is  with  respect  to  P coordinates, 
and  requires  differentiation  with  respect  to  the  implicit  dependence  of 
F on  T|  through  the  dependence  of  F on  t(^),  and  also  with  respect 
to  the  expl  icit  dependence  of  F on  T).  The  partial  s ur face  d i verge nee 
of  F is  defined  by 


F=f-VF'P+m*VF*m  , 

P P 


where  the  partial  gradient  operator  is  with  respect  to  the  explicit 

dependence  of  F on  P only,  with  the  t (T))  component  of  F held 
fixed.  The  total  surface  divergence  of  F defined  in  (2.4)  can  then 
be  written 


7s  • F • V ■ V * %S  • f 


The  surface  gradient  of  t is  defined  by 


= Vt  - (n  • Vt)  n , 


1 


6. 


and  i s 


_ m sin  8 

V * JTT^iTT 


Total  surface  divergence  of  F is  therefore 


. f = 3F  . HI  sin  G . - . 

S c(l-Mr)  PS  • y 

This  identity  plays  a fundamental  role  in  what  follows. 

There  is  a theorem  in  differential  geometry  that  says  that  the 
integral  over  a closed  surface  S of  the  total  surface  divergence  of  a 
vector  F everywhere  tangent  to  S is  zero  (Weatherburn,  Reference  5a) : 


V$  * F dS  = 0 , F tangent  to  closed  S 


(2.6) 


If  F is  tangent  to  S , integration  of  (2.5)  over  S gives 


' *F 

J Trr 


m s i nG 
cO-M  ) 


dS  = - 


'ps 


F ds 


(2.7) 


The  theorem  must  be  understood  with  the  restriction  that  F is  not 
too  singular  on  S ; this  consideration  plays  a guiding  role  in  what 
fol lows . 


The  integrand  on  the  left-hand  side  of  (2.7)  has  a Doppler 
factor,  while  that  on  the  right-hand  side  is  free  of  this  factor.  If 

F is  a reasonably  well  behaved  function  on  S , then  V • F can  never 

PS 

produce  singular  behavior  on  S because  the  t (fj)  component  of  F is 
held  fixed  in  this  operation.  However,  what  is  required  is  not  the 
integrand  that  appears  on  the  left  of  (2.7),  but  rather  an  integrand 
of  the  form  F cos0/(l-Mr)  that  appears  in  the  fundamental  radiation 
inteqral  of  (2.2).  The  identity  (2.7)  must  be  further  defined  so  that 
the  basic  radiation  integral  of  (2.2)  emerges  in  terms  of  less  singular 
integrals . 
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3.  SIMPLE  IDENTITIES 


Equation  (P.5)  can  be  rewritten  identically  as 

vs  ' "F  = cO-mT  ^ ‘ "F  sin0)  - cdT)  sin0)  + VPS*  F • {3J) 

The  t - derivative  of  m sin0  is 


±,~ 


(m  s inO)  = ? - (n  * ? ) n . 


The  super  dot  means  */9t  at  constant  T)  and  r = - <u  x r . For  any 
function  F , scalar  or  vector,  we  have  the  identity 

^f  = $£ 

St  ^t  _ 1-M 


(3  .2) 


(3.3) 


Use  of  (3.2)  and  (3.3)  in  (3  .1  ) gives 

VF  = ^ (™  • F sin6)  • c(V-Mr)  ’ f n - (n  • r )n]  + Vps  ‘ F . (3.4) 

This  identity  is  made  useful  in  radiation  integrals  by  defining  the 
vector  F by 

F - X(t(1M)  F (T(1j)f1j)  , 

where  now  X and  F are  an  arbitrary  vector  and  scalar  defined  on  S . 

We  write  X in  component  form  as  £ = arf  + £m+Yn.  The  n 
component  of  X cancels  throughout  equation  (3.4).  We  therefore  redefine 
X by  X = a x n , where  5 is  arbitrary.  Equation  (3.4)  takes  the  form 

Vs'Orx  n F = - ^-(F  a x m • n sin9)-  (c  (i-M  )°'  x n * ^ + VPS  * a xn  F • (3  .5) 

The  vector  a can  be  written  in  component  form  as  a = AT  + Bm  , where 
A and  B are  arbitrary  scalars  defined  on  S.  Equation  (3.5)  becomes 


sc m 
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- v _ A /FA  sln0*  . F 
7S  (er-Am)F  = - — ( ) + cTi-M  ') 


+ 7ps  • (Bf  - Am)F 


(Auj  • f cos0  - Boj  " N) 


(3.6) 


The  unit  vector  N = T x r . Some  simplification  results  if  we  write 
A and  B as  A=a,  B=-b  cos©  . The  identity  (3.6)  is  then 

(br  cos0  + am)  F = ( ) - ^jTj^-)  ^ * " + b au  * N) 

+ (8  cos0 f + am)  F . 


(3.7) 


This  equation  is  actually  two  equations  determined  by  assuming  a / 0, 
b = 0,  and  then  a = 0,  b 4 0.  We  require  the  combined  form  (3.7). 

The  identity  (3.7)  is  used  to  solve  for  a term  of  the  form  F cos0/(l -M  )"' . 
To  do  so,  define  a new  scalar  function  G by 

G = F(au>  • f + b£  * N)  . 


This  factor  appears  in  the  second  term  on  the  right-hand  side  of  (3.7). 
In  terms  of  arbitrary  scalars  a,b,  and  G,  equation  (3.7)  is 

^ . T (bfcos0  + am)  G ~|  _ 3 |~  Ga  s i nB "]  G cosB 

S aJ  * f + bu>  * N " c (aui  * f + boi  ’ N)  C ^‘Mr^ 


. r (bF  cosQ  + am)G 

ps  L - 


aou 


f + bi  • N) 


(3.8) 


Integration  of  (3.8)  over  S and  application  of  the  integral  theorem 
(2.6)  gives  the  identity 


l*G  cos0  <3  T Ga  sin© 

J "7T— m”T  dS  = 37  J 

' r aou  r + bo 


dS 


bai  ’ N 


+ c 


f . I-  (bf  cos0  + am)G 

J PS  L 

am  f + bm  N 


dS 


(3  .9) 


r 
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The  scalars  a and  b In  equation  (3.9)  are  still  arbitrary  and 
at  our  disposal.  It  is  evident  that  a variety  of  identities  can  be  generated 
by  various  choices  of  these  scalar  functions.  One  identity  that  is 
particularly  useful  is  based  on  setting 

a = ou  ' f = (Up  , b = uj"N=u)N  . 


The  denominators  in  (3.9)  are  then 

auip  + bu>N  = (ii  • f)2  + (i  * N)2  . (3.10) 

The  unit  vectors  f , r and  N = P x r are  an  orthogonal  system 
defined  at  each  point  of  S . The  angular  velocity  ou  can  therefore 
be  written  as 

£«(«*  f)  f + (®  *?)?+(£  * N)N  . 


Then  (3.10)  Is 


aojp 


* b*N- 


ou 


(£  • r)‘ 


Let  us  call  this  quantity 


n 


2 

to  - (ou 


r)2 


(3.11) 


The  identity  (3.9)  becomes 

I G(!-h1  ds  “ I n"2  G "V  slne  dS^JvPs[n‘2(-/cose  + V;)G>s  (3  A2) 

' ~r 

Now  1 - derivatives  of  ? with  t held  constant  are  all  0(l/r),  and 

therefore  so  are  1)  - derivatives  of  fl  with  t constant.  Furthermore, 

_2 

the  t derivative  of  fi  Is  Identically  zero.  Hence  Cl  Is  either 
asymptotically  or  identically  zero  with  respect  to  all  differential 
operations  In  equation  (3.12),  and  this  factor  may  be  treated  as  a 
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constant  factor  In  that  equation.  Thus  we  finally  have  the  identity 

J dS  = n”?  St  I G “>  x dS  + cn"2  j f cosO+aipin)  G^JdS  . 

(3.13) 

The  relation  o>p  sin®  = ou  x r * n has  been  used  for  the  second  integral 
of  (3.13).  The  left-hand  side  of  this  equation  has  the  form  required 
in  the  radiation  Integral  (2.2).  However,  (3.13)  can  be  further 
simplified  in  a useful  way.  The  vector  in  the  last  integrand  of  (3.13) 
can  be  wr i tten  as 

CD^5  cos©  + cupm  = • id  x n + (J  ■ r)  r x n . (3.1*0 

This  identity  is  easily  verified  by  taking  the  scalar  product  of  both 

sides  successively  with  f and  m.  The  last  integrand  of  (3.13)  is 

7PS  * [ (u,Nr  cos6  + <Y")G  J = cos6  + «V>  ^rsG 

+ G Vps  * (mj  cos0  + ciym) 

The  last  divergence  in  this  equation  is,  using  (3.1*0, 

7PS  * (“1/  cose  + tY"')  * - 7PS  * (Sx  n)  ♦ Vps*  f*(J  • f)rx  IT].  (305) 

The  angular  velocity  a>  is  constant  in  (3.15)*  and  the  T]  - derivatives 
of  the  factor  ou  * ? are  0(l/r)  because  l)  - derivatives  with  T 
constant  of  r are  0(l/r).  Hence  equation  (3.15)  Is 

7ps  • («,,/  cos0  + cop^)  « S • 7^  x n . (»  • p)f  • 7p$  x n . (3  .16) 

But  the  surface  curl  of  the  unit  normal  surface  vector  n is  zero; 

this  Is  a fundamental  relation  of  differential  geometry  shown,  for  example, 


In  Reference  5».  Hence  equation  (3.16)  Is 

7PS  ’ COS°  + (V’)  * 0 • 

and  the  third  Integrand  of  (3.13)  is 

Vp$  (ttJ^T  cose  + uiprn)G  ■ ■ i x n ' 7p^G  + (u>  * r)rxn  * 7p$G  . 

The  fundamental  identity  (3.13)  is  therefore 

J T-ir* ds  E n"2  St  J G ® x ? ' " dS 

r 

- c n"2  J J x n • Vp$G  dS  + (w  • r)CT2  Jr  x n * 7p$  G dS  (3.17) 


We  now  identify  the  function  G in  (3.17)  with  the  source  strength 
F/r  in  the  radiation  integral  of  equation  (2.2).  The  contribution  to 
far  field  pressure  from  sources  distributed  on  S is 


>2  r c - * - ds  n-2  ^ r 

, Fwx  r . n-_-n  -J 


V F^ 
PS  r 


+ CT2(w  * r)  J f 


V F — 
PS  r 


If  the  observer  is  in  the  plane  of  rotation,  ui  * r = 0 and  the  last 
Integral  In  (3.18)  Is  zero.  If  the  surface  S Is  hovering  (translational 
velocity  v - 0 ),  the  source  strength  F is  independent  of  t in 
coordinates  fixed  to  S . The  integrand  of  the  second  integral  in 
(3.18)  Is  therefore  Independent  of  t and  hence  Independent  of  t for 
hover,  and  the  time  derivative  of  this  Integral  is  zero.  Therefore  for 
hover  and  observer  In  the  plane  of  rotation 

, 0) ^ r - - * ~ dS  /,  . n\ 

4rrp  - — £•  F ui  x r • n — . (3.19) 

c ^ J r 


l 
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The  inteqral  depends  on  observer  position  and  time  throuqh  the 
dependence  of  r on  T(“,x,t).  Equation  (3.19)  appears  to  be  a remarkably 
simple  expression  for  this  important  case. 


4.  GEODESIC  PROPERTIES  OF  THE  ROTATING  SURFACE 
The  last  inteqral  in  equation  (3.18)  is 


r * - . _ c dS  r • Q p.  „ r dS 

. r x " 7p$F  T-  J sm0  . Vp$F  — . 


This  result  follows  from  the  definition  of  the  unit  vector  f in 
equation  (2.3).  For  any  scalar  F , 


P . _ *F  f 'm  sinP  . p . _ 

r V = V c7i-HrV  + r 7PSF  • 


Since  f • m = 0,  we  find  f • 7^F  = f * 7p<_F.  Equation  (4.1)  is 


therefore 


J f x n • V^F  ^ = J sinP  f * V$F  ^ . 


For  the  second  inteqrand  in  (4.3),  we  have 


sinP  f • 7$F  s'  7S  • (fF  sinP)  - F 7$  • (fsinP)  . 


The  factor  r can  always  be  treated  as  constant  in  derivatives 
with  respect  to  T - coordinates  and  t . Therefore  (4.3)  is 


r x n 


The  second  inteqral  is  zero  by  the  surface  inteqral  theorem  (2.6) 
In  the  third  inteqral  of  (4.4), 


7^  * T sinP  * sinP  7^  * f + cosP  f * 7 0 


Now  the  unit  vectors  f and  m form  an  orthogonal  net  on  S.  The 


total  surface  diverqence  of  f is  1 • f = V,  where  V is  the 


But  see  the  note  added  in  proof,  page  21. 


(4.1) 


(4  .2) 


(4.3) 


(4.4) 


(4.5) 


I 
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qeodesic  curvature  of  the  coordinate  curve  on  S with  m as  unit 
tanqent  vector.  The  directional  derivative  f ' Vr9  in  (4.5)  is 
f • 7^9  = f * Vp^P,  an(l  this  equality  results  from  (4.2)  and  f • m = 0. 
Since  n • r = cos9, 

- sinP  f • 7p?9  = T * Vps  n 1 r = n • r , (4.6) 

where  is  differentiation  alonq  a F - curve  with  t constant. 

The  last  term  in  this  equation  is 


n ' r = r • = - r • (h^T  + T m ) , (4.7) 

where  T is  the  torsion  of  a qeodesic  curve  tanqent  to  a f - curve, 

and  h is  normal  curvature  of  an  m - curve.  The  last  equality  in 
n 

(4.7)  follows  from  differential  properties  of  the  unit  normal  vector  n 
as  discussed  in  Chapter  I of  Volume  II  of  Wea therbur n 1 s texts.  Orthogonality 
of  r and  f reduces  (4.7)  to 

— n • r = - T sinP  . 

Use  of  this  relation  in  (4.6)  qives 


f • ’os  “ ‘ T • 


and  hence  equation  (4.5)  is 


7^  • I*  s in'5  = V s in0  + T cos  9 . 


The  left-hand  side  of  (4.4)  can  therefore  be  written 

| r x n 1 7p^F  = - J (v  sin9  + T cos9)F  ~ 


(4.8) 


Surface  derivatives  of  F in  the  first  inteqral  of  (4.8)  have  been 
eliminated.  This  result  is  accomplished  at  the  expense  of  introducing 
the  relatively  unfamiliar  qeodesic  properties  of  S . 
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Use  of  equation  (4.8)  in  (3.18)  qives  an  alternative  expression  for 
acoustic  pressure  due  to  surface  integrals: 


4rrp  = 


Q~?~  > 

C ^t2 


r - - . - ds  n-P  > n r dS 

F (D  x r • n — - n — (1)  x n V_F  — 
r *t  J PS  r 


- 0 2 (u)  ‘ r)  ^ J (V  sinP  + T cosP)F  ^ 

O 

Throughout  the  preceeding  formulas,  the  factor  0 is 

zero  when  the  retarded  radiation  vector  r is  parallel  to  the  anqular 

velocity  id  . This  occurs  when  ? is  perpendicular  to  the  translational 

velocity  V and  the  anqular  velocity  id  x r . Retarded  times  t over 
2 

which  u is  zero  or  small  therefore  correspond  to  values  of  t over 
which  the  radiation  Mach  number  M is  zero  or  small.  This  situation 


(4.9) 


r 

occurs  when  the  lateral  distance  of  the  observer  from  the  center  of 
rotation  is  held  fixed  and  the  vertical  distance  becomes  infinite  in  the 
limit  r -*  ®.  This  observer  position  corresponds  to  the  rotary  surface 
beinq  straight  overhead  (or  below)  at  some  retarded  time.  In  this 
case,  acoustic  pressure  can  be  determined  by  radiation  integrals  in 
which  M^  is  zero  or  small,  and  the  generalized  Lowson  formulas  are 
well  suited  to  that  kind  of  calculation.  Equations  (3.18)  and  (4.9)  are 
useful  for  those  observer  positions  at  which  the  amplitude  of  radiated 
sound  may  be  large.  These  positions  may  be  in  the  plane  of  rotation, 
or  on  any  line  through  the  center  of  rotation  that  is  inclined  at  an 
angle  of  less  than  90°  to  the  plane  of  rotation. 


5 . SHOCKS 

1.  The  hovering  rotor. 

The  shock  wave  on  a hovering  surface  S rotating  at  constant 
angular  velocity  Is  stationary  relative  to  S.  The  geometry  of  the  shock 
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is  independent  of  time,  and  the  velocity  relative  to  an  Inertial  frame 
of  points  on  the  shock  is  solenoidal  . The  hoverinq  shock  wave  may  therefore 
be  considered  as  a riqid  surface.  This  description  of  the  shock  surface 
is  no  longer  strictly  true  if  the  body  surface  is  rotatinq  and  3lso 
ad vane i nq . 

Let  S'  be  the  shock  surface  on  a hoverinq  rotor.  Let  n ) be 
the  unit  vector  normal  to  S*  that  points  in  the  direction  upstream  of 
S,t  and  let  0 be  the  anqle  between  n and  the  unit  radiation  vector 
? (r ) . The  contribution  to  acoustic  pressure  from  the  shock  is 


^ f dS ' 

krrep  = - J ATrr  cose 


Integration  is  over  the  shock  surface  and  the  integrand  is  evaluated 

at  retarded  time.  The  source  strength  AT  is  the  jump  in  the  resolved 

Lighthill  stress  tensor  T across  the  shock.  If  the  superscript  (2) 

denotes  the  upstream  side  of  S'  and  (1)  the  downstream  side,  then 

AT  = T - T * , and  T = r .? . T.  . . 
rr  rr  rr  rr  i j ij 

Equation  (5.1)  is  the  equivalent,  in  surface  fixed  coordinates, 
of  a similar  contribution  found  by  Farassat  in  his  coordinate  representation, 
However,  he  finds  an  additional  term  that  involves  the  derivative  in  the 
radiation  direction  of  ATr(_.  That  term  appears  to  be  in  error  and  is 
omitted  here. 

The  Doppler  factor  in  (5.1)  can  be  suppressed  by  transforming  the 
integral  along  the  lines  used  in  sections  2 and  3>  only  now  the  shock 
surface  S*  is  not  closed  and  a more  general  version  of  the  surface 
divergence  theorem  of  equation  (2.6)  is  required. 

Let  C be  the  closed  curve  that  bounds  the  shock  surface  S*. 

Let  / be  a unit  vector  normal  to  C and  tangent  to  S'  and  pointing 
away  from  S7.  Let  u be  a unit  vector  tangent  to  C so  that  5, u and 
n from  a right  handed  system  of  unit  vectors  on  C . If  F(t(1)),^)  is 
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any  vector  everywhere  tanqent  to  S,t  then  (Weatherbur n,  Vol  . I) 

j*7s*F  dS/  = Jf*2  ds  . (5.2) 

We  aqain  introduce  unit  vectors  F and  m tanqent  to  S'  as  defined 
by  equations  (2.3).  with  total  surface  divergence  r’^,F  defined  by  equation 
(2.4).  The  contour  integral  on  the  riqht-hand  side  of  (5.2)  is  taken  around 
C in  the  positive  direction  determined  by  the  tanqent  vector  u ; ds  is 
arc  lenqth  alonq  C . 

Trans  format  ion  of  ^5.1)  now  follows  the  same  steps  used  to  reduce 
integrals  over  closed  surfaces,  developed  in  sections  2 and  3.  The  basic 
surface  divergence  identity  for  an  arbitrary  scalar  G is  equation  (3.8), 
which  can  be  written  as 


c o'pvr-r (u)^ 


-?  3 


s cosP  + (Upm)G]  = 0 ' ^-(G  .Up  sinP) 


" + c n'2vps’t  (wNr  cosfi  + , 

r 

2 2 - * 2 

with  Q = a)  - (<!>•?)  . The  vector  in  parantheses  in  (3.8/)  is  again 


(3.8' 


wj  cosO  + = - mxn  + (wr)rxn  . (3.I47) 

The  partial  surface  divergence  in  (3.87)  is  aqain 

7ps-r  (u^r  CO sP  + OY.m)G]  = (wj*  cosO  + iu^_m)  Tp5G 

= - Jxn*VpsG  + (i‘f  )fxn-VpsG  , (5.3) 

using  (3.l4,)in  the  last  equality.  Now  substitute  (3.14')  and  (5.3) 

in  (3.8/),  integrate  over  S7  and  apply  the  open  surface  integral  theorem, 

equation  (5.2).  Identify  the  arbitrary  function  G in  (3.87)  with  the 

shock  source  strength  VT^/r  in  the  radiation  integral  of  equation  (5.1). 

Note  that  the  shock  surface  5 7 generally  intersects  the  body  surface 

S in  a curve  Cj  that  is  part  of  the  shock  boundary  C.  The  part  of 

the  boundary  C that  does  not  lie  along  the  body  surface  is  just  the 

sonic  line  of  the  shock.  On  this  part  of  C , AT  * 0.  The  closed 

r r 

countour  integral  in  (5.2)  therefore  reduces  to  an  integral  along  the 


SBC®*' 
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intersection  of  the  shock  S*  with  the  body  S»  and  this  is  a line 
inteqral  alonq  C(.  The  contribution  to  acoustic  pressure  from  a 
shock  is  therefore 


4tt  0 p = 


*t 


2 J 


f id*  r • n 


AT 


dS' 


rr  cr 


A P - - 

M J 


VDC  AT 

PS  rr 


dS' 


(5  .4) 


For  a hovering  rotor,  the  second  and  fourth  integrals  on  the  right-hand 
side  of  (5.4)  are  independent  of  t because  their  integrands  are  in- 
dependent of  retarded  time  t . In  this  case,  time  derivatives  of  these 
inteqrals  are  zero.  If  the  observer  is  in  the  plane  of  rotation,  u>-r  = 0. 
Therefore,  for  hover  and  an  in-plane  observer,  the  shock  contribution  to 
radiative  pressure  is 


-?  v r — „ — ds 

4rrp  = (D  — — ouxr'n  AT  

M2  J rr  cr 

This  formula  is  analogous  to  equation  (3.19)  that  gives 
to  acoustic  pressure  from  a closed  hoverinq  surface  and 


(5.5) 

the  contribution 
an  in-plane  observer. 


2.  Advancing  rotor. 

If  the  forward  velocity  of  the  rotor  is  not  zero,  the  shock  surface 
is  not  necessarily  stationary  relative  to  the  body  surface.  In  addition, 
the  shape  of  shock  surface  may  chanqe  as  it  grows  and  collapses  with  the 
passage  of  the  blade  tip  throuqh  the  advancing  azimuthal  region.  The 
deformation  of  the  shock  means  that  the  velocity  of  points  on  it, 
relative  to  an  inertial  frame,  is  no  longer  solenoidal  . Throughout 
But  see  the  note  added  In  proof,  page  21  . 
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Reference  3 it  was  assumed  that  both  the  body  and  shock  surfaces  are  rigid, 
and  hence  that  their  velocities  have  zero  divergence.  It  may  be  that  a 
changing  shock  geometry  produces  only  minor  variations  in  far  field 
pressure,  but  the  assumption  of  a rigid  shock  geometry  is  probably  not 
quite  correct,  and  possible  errors  introduced  by  this  assumption  are 
somewhat  obscure.  The  results  of  Reference  3 could  be  altered  to  account 
for  unsteady  shock  geometry;  the  changed  analysis  would  begin  with  equation 
(5)  of  that  paper  and  would  apply  only  to  contributions  from  the  jump  in 
the  Lighthill  stress  tensor  across  the  shock.  It  is  also  assumed  in 
Reference  3 that  the  sign  of  the  local  first  curvature  of  the  shock  does 
not  change  from  point  to  point. 

Expressions  for  the  contribution  to  acoustic  pressure  from  an 
advancing  shock  are  greatly  simplified  if  the  shock  can  be  considered  as 
rigid,  and  this  assumpt ion  will  be  made  here.  The  boundary  C of  the 
shock  surface  s'  must,  of  course,  be  considered  as  a function  of  time. 

We  shall  also  assume  that  the  acceleration  of  the  shock  surface  relative 
to  the  body  is  small,  or  that  variations  in  this  relative  velocity  occur 
over  such  a small  time  inteval  that  they  may  be  neglected.  This  assumption 
means  that  the  T-derivative  of  the  unit  radiation  vector  ? is  still 
? * - <dx?  and  the  fac  or  C defined  by  equation  (3.11)  can  be  regarded 
as  independent  of  t and  ^ coordinates.  Assumptions  about  shock 
geometry  and  acceleration  can  be  removed  at  the  expense  of  added 
complication.  The  extent  to  which  these  assumptions  are  reasonable  may 
be  clarified  when  we  have  detailed  calculations  that  show  shock  behavior 
on  advancing  blades. 

With  the  above  assumptions,  equation  (5.4)  gives  the  contribution 
to  acoustic  pressure  from  the  unsteady  advancing  shock  wave.  The  second 
integral  in  equation  (5.4)  can  be  rewritten  using  equation  (4.8).  Surface 
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derivatives  of  AT  are  thereby  eliminated  from  the  inteqral,  and 
rr 

far  field  pressure  will  depend  on  local  qeodesic  properties  of  the  shock 
s ur  face . 


5.  SUMMARY  AND  CONCLUSIONS 

Fquation  (3.18)  gives  the  contribution  to  acoustic  pressure  from 
a distribution  of  sources  of  strength  F on  the  rotating  and  translating 
surface;  the  corresponding  contribution  from  shock  waves  is  equation  (5.4). 
If  the  surface  is  hovering  and  the  observer  is  in  the  plane  of  rotation, 
these  formulas  simplify  to  their  respective  counterpar ts , equations  (3.19) 
and  (5.5).  Each  radiation  inteqral  is  free  of  singularities  as  long  as 
the  surface  advancing  tip  Mach  number  is  not  supersonic.  The  integrals 
are  a function  of  aerodynamic  data  on  the  blade  and  shock  surfaces,  and 
this  Information  is  assumed  to  be  known  from  near  field  calculations.  The 
near  field  numerical  calculations  are  always  carried  out  in  a coordinate 
system  fixed  to  the  accelerating  surface.  The  required  data  can  therefore 
be  inserted  in  radiation  integrals  with  a minimum  of  transformation. 

The  source  strength  F in  equations  (3.18)  and  (3.19)  is  qenerally 
a linear  combination  of  the  Lighthill  stress  tensor  (resolved  twice  in  the 
radiation  direction),  the  pressure  distribution  on  the  blade  surface  and 
also  the  surface  radiation  Mach  number.  The  complete  far  field  pressure 
is  composed  of  body  surface  and  shock  contributions  according  to  equations 
(3.18)  and  (5.4);  in  addition,  it  depends  on  the  derivative  in  the  radiation 
direction  of  the  Lighthill  stress  tensor,  which  must  be  integrated  over  the 
volume  exterior  to  the  body  and  shock  surfaces.  This  residual  volume 
integral  is  not  discussed  here.  Its  inteqrand  is  a second  order  term,  and 

the  integral  may  be  negligible. 
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NOTE  ADOED  IN  PROOF 
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The  contribution  to  acoustic  pressure  from  sources  distributed 
on  a hovering  surface  is  given  by  eguation  (3.19);  the  contribution 
from  the  shock  on  the  hovering  surface  is  equation  (5.5).  Both  of  these 
results  are  based  on  the  assumption  that  the  source  strengths  are  independ- 
ent of  time,  and  therefore  independent  of  retarded  time,  in  coordinates 
fixed  to  the  surface.  This  assumption  is  false,  or  rather  it  is  half 
true,  for  both  the  body  surface  and  the  shock.  The  body  surface  source 
strenqth  F in  equation  (3.19)  generally  depends  on  surface  pressure 
and  density,  and  both  of  these  quantities  are  independent  of  retarded 
time  in  7)  coordinates.  However,  E also  depends  on  the  component  of 
the  disturbance  velocity  in  the  radiation  direction,  and  this  quantity 
is  a function  of  retarded  time  r through  the  dependence  of  the  radia- 
tion vector  r on  t . Thus  F has  components  for  which  equation  (3.19) 
is  correct,  and  also  components  for  which  the  full  equation  (3.18)  is 
required.  The  same  must  be  said  of  the  shock  source  strength  AT  : 
it  has  components  that  depend  on  the  jumps  in  pressure  and  density  across 
the  shock,  and  both  of  these  quantities  are  independent  of  T in  71 
coordinates.  Equation  (5.5)  is  correct  for  the  contribution  from  these 
discontinuities.  However,  ATrr  also  has  components  dependent  on  the 
disturbance  velocity  resolved  in  the  radiation  direction  r , and  these 
contributions  to  far  field  pressure  must  be  expressed  by  the  full 


equat ion  (5  .4)  . 


